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Abstract 

An analytic theory is presented for the width selection of Saffman- Taylor fingers in the 
presence of thin film effect. In the limit of small capillary number Ca and small gap to 
width ratio e , such that e << Ca « 1 , it is found that fingers with relative width 
A < \ are possible such that 7* > where tlie P°si tive constant k depends 

on the branch of solution and equals 2.776 for the first branch. A fully nonlinear analysis 
is necessary in this problem even to obtain the correct scaling law. It is also shown how in 
principle, the selection rule for arbitrary Ca can be obtained. 
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I. Introduction 


In recent years, considerable theoretical work, both numerical (see McLean & Saffman 
(1981), Romero (1982), Vanden-Broeck (1983), Kessler & Levine (1985), Tanveer (1986) 
and Tanveer (1987a)), and analytical (see, for example, Combescot et al (1986), Shraiman 
(1986), Hong & Langer (1986), Tanveer (1987b), Combescot et al (1988), Dorsey & Martin 
(1987), Comescot &: Dombre (1988), Hong & Langer (1988) and Tanveer (1989a)) has been 
focussed on the singular effect of surface tension that breaks the continuum of Saffman- 
Taylor (Saffman & Taylor (1958), Taylor & Saffman (1959)) exact solutions for steadily 
propagating fingers and bubbles through a viscous fluid in a Hele-Shaw cell under the 
assumption that the thin film of the more viscous fluid in the narrow gap direction plays 
no significant role and that there is no variation of curvature in the transverse (narrow gap) 
direction. The above simplifying assumptions leads to a set of relatively simple boundary 
conditions on the flow variables at the interface and will be referred to as the McLean- 
Saffman (MS) boundary conditions henceforth. In this simplification, the mathematical 
equations for a finger contain a single non-dimensional parameter Q = , where 

b is the gap width and 2 a the cell width of the Hele-Shaw cell; T , fj, and U are the 
surface tension, viscosity of the displaced fluid and the finger velocity respectively. The 
analytical and numerical evidence to date suggests that with the MS boundary conditions, 
for Q over some range, a discrete set of steady solutions are possible, each characterized 
by a different velocity. In the limit of zero Q , the finger or bubble velocity U — ► 2 V for 
all branches of solution, where V is the velocity of displaced viscous fluid at infinity. For 
the finger, with the MS boundary conditions, the continuity of fluid flow implies that its 
relative width A = V/U , from which follows that as Q — ♦ 0 , A — > ^ . It is also 
known that the selection of the finger (Combescot et al (1986), Shraiman (1986) and Hong 
h Langer (1986)) or bubble velocity (see Combescot & Dombre (1988), Hong & Family 
(1988) and Tanveer (1989a)) as well as its symmetry (see Tanveer (1987b) and Combescot 
& Dombre (1988)) about the channel centerline for small Q is caused by transcendentally 


2 



small corrections in Q that are beyond all orders of an asymptotic series in powers of Q . 


It has also been shown that such terms play a crucial role in the linear stability (see 
Tanveer (1987c), Bensimon et al (1987), Tanveer (1989a)) as well (for the MS boundary 
conditions). When Q = 0 , there is a continuum (see Tanveer (1987c) , Tanveer & 

Saffman (1987)) of unstable modes for the finger or the bubble suggesting that the time 
dependent problem with zero surface tension is ill-posed. However, if Q is small but 
non zero, there are in general no eigenvalues close to the eigenvalues of the zero Q lin- 
ear stability operator. In the small Q limit, it has been shown that this is due to the 
effect of transcendentally small terms in Q . Numerical (see Kessler & Levine (1986), 
Bensimon (1986) and Tanveer & Saffman (1987)) and analytical (see Tanveer (1987c) and 
Tanveer (1989a)) calculations show that only one branch of finger solutions is stable, while 
others are unstable for arbitrary Q . However, numerical simulation by Degregoria & 
Schwartz(1986) of the time dependent problem based on the MS boundary conditions 
show instability at sufficiently small values of the surface tension parameter, despite the 
linear stability of these solutions. It is now understood (Bensimon (1986) that the insta- 
bility observed in the time dependent calculations is a manifestation of all the branches 
of solution tending to the same A = ^ Saffman-Taylor solution. Since all but one is 

unstable, the threshold amplitude of instability must tend to zero as surface tension tends 
to zero. Bensimon (1986) further conjectures that the threshold amplitude of instability 
tends to zero exponentially with surface tension and his numerical results tends to support 
this though accurate verification of the claim does not appear possible by direct numerical 
calculations. It has been suggested (Kessler &: Levine, 1988) that the different branches of 
solution differ by an exponentially small amount as Q tends to zero and this would there- 
fore explain the exponential dependence of threshold amplitude of instability. However, 
we disagree with this explanation since A and the shapes on different solution branches 
differs by terms of order C/ 2/3 for small Q , which is not consistent with exponential 
dependence of the threshold level of instability, if indeed it is true as Bensimon suggests. 
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Aside from the point about the threshold level of instability raised in the last para- 
graph, the behavior of solutions to mathematical equations based on the McLean- Saffman 
boundary conditions appear to be quite well understood. However, experiments by Saffman 
& Taylor (1958), Tabeling et al (1987) and Kopf-Sill &: Homsy (1987) with the finger show 
quantitative discrepancies with the theoretical finger width predictions of the MS theory. 
The bubble experiments (Maxworthy (1986), Kopf-Sill &: Homsy (1988)) also show signifi- 
cant disagreements. Detailed experiments with the finger by Tabeling et al (1987) suggest 
that there isn’t a single control parameter Q as would be the case if the MS theory were 
accurate. Experiments of Tabeling et al (1987) also show that to a significant extent, the 
gap to width ratio e = ^ and capillary number Ca = ^ are two independent control 
parameters. Equivalently, one can think of Q ( = ) and Co as two independent 

parameters. This suggests that the thin film effects play a significant role in the fingering 
problem that have to be accounted for to explain discrepancies between the MS theory 
and experiment for the the finger width prediction. There is also a significant discrepancy 
on the onset of instability for small Q . Experiment (Tabeling et al, 1987) shows that the 
finger is unstable for small Q . However, there does not appear to be a sharp critical value 
of Q below which the finger is unstable; the instability point appearing to depend on the 
noise level in the experiment. This suggests perhaps a nonlinear instability mechanism 
(at least under some experimental conditions) as suggested by Bensimon (1986) for the 
MS boundary conditions. However, the threshold noise level that destabilizes the finger 
depends both on Q and Ca suggesting that the Bensimon mechanism for non-linear 
stability needs significant modifications if it is to be valid for boundary conditions that 
include thin film effects. 

Saffman (1982) discusses the general form of the boundary conditions necessary to 
include the thin film effects. Park &: Homsy (1985) and Reinelt (1987a) further develop 
details of these conditions that incorporate thin film effects into gap averaged 2-D bound- 
ary conditions for the steady finger, which we call Saffman-Park-Homsy-Reinelt (SPHR) 
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conditions. In addition to the pressure drop across the interface due to the lateral curva- 
ture of the interface that is included in the MS boundary conditions, the SPHR conditions 
include an additional pressure drop term due to variation of transverse curvature at dif- 
ferent points on the interface in the lateral plane. It also includes terms in the kinematic 
boundary conditions that account for flow into the thin film region. The SPHR conditions, 
as we shall see in this paper, are quite complicated since they involve knowledge of func- 
tions that have to be determined by solving an associated 3-D problem in the transverse 
plane. In general, the solution to this associated three dimensional problem has to be com- 
puted numerically. For small Ca , using Bretherton (1961) results, Park fc Homsy (1985) 
and Reinelt (1987a) have found asymptotic expressions for these functions. Numerical 
calculations by Reinelt (1987b) and Sarker & Jasnow (1987) show that the discrepancy 
between theory and experiment on the finger width dependence on the control parameters 
is greatly reduced when the SPHR boundary conditions are used. Equivalent boundary 
conditions for the bubble or for the time dependent problem are yet to be deduced. 

The numerical work incorporating thin film effects cited in the last paragraph is impor- 
tant; however it leaves some other important questions unresolved. First is the existence 
of other branches of solutions. Reinelt (private communication) has found some other 
branches over some parameter ranges. For reasons stated earlier, it is important to know 
if all these steady state solutions tend to the same width as Q — > 0 , or if it is possible 
for the limiting finger width to be different for different branches. Also, the precise de- 
pendence of A on Q and Ca in the limit of small Q is difficult to conclude from the 
numerical solutions for we now have two control parameters. To compound our problem, 
the problem is numerically ill posed as Q — * 0 for any Ca as a continuum of solution 
with arbitrary width exists for Q = 0 and, as seen in this paper, transcendent ally small 
terms in Q have to be resolved for determination of A . The need of an analytic theory 
to resolve this limit is therefore obvious. 

The purpose of this paper is to incorporate the thin film effect into an analytic theory 
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which is valid for small Q . Further, we restrict ourselves primarily to the case of small Ca , 
where simplifications of the SPHR boundary conditions are possible by using Bretherton’s 
formulae, though in section 7, we show how the analytic theory can be extended for any 
Ca by some subsidiary numerical calculations. Also, our primary concern will be the 
possibility of solution for A < ^ , though we show that selection of finger width with 

A > ^ is possible for a few restricted cases in the two-parameter space. We also present 
conditions under which the selection results of the theory based on MS boundary conditions 
hold. However, we do not present an exhaustive study of selection for all possible relative 
orderings of the small parameters Ca and Q . 

For A < ^ > we conclude that it is possible for solutions to exist for Q « Ca << 
1 with A 1-V^ = k , where the positive constant k depends on the branch of 

solution and equals 3.061 for first 1st branch. In terms of control parameters e and Ca , 
this implies that for e << Ca <<1, - = — * / a , where the positive constant 

k depends on the branch of solution and equals 2.776 for the first branch. The lower limit 
on Ca for the validity of the above relation is sufficient to ensure that the solutions found 
satisfy the assumptions made in the analysis. Whether the relation hold over a larger range 
of the parameter space is an open question. We notice that when Q l/2 « Ca « 1 , 
A is close to zero. On the otherhand if 1 >> £ 1/2 >> Ca » Q , A is close to a 
half. If Q ^ 2 = O(Ca) , for sufficiently small values of the two parameters Ca and Q , 
one can obtain fingers with width in the open interval (0, 5) . In this range of Ca and 
G , the A on different branches do not tend to the same value as Ca , Q — ► 0 . 

For A > 5 , we conclude that the finger solutions calculated on the basis of the MS 
theory for which A ~ ^ + constant Q 2/S persists when Ca << Q 7/,a << 1 . This is 
rather unexpected since on inspection of the SPHR boundary conditions and their limiting 
form using Bretherton’s formulae, it would appear that thin film or transverse curvature 
term or both (call them the 3-D effects) are significantly bigger than the lateral curvature 
term that is included in the MS theory when Ca » Q 3 . Thus for Q z « Ca « g 7 ' 3 , 
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we have a situation where to the leading order, the solutions based on the MS theory gives 
the correct finger shape and the width scaling with Q even when the 3-D terms in the 
SPHR boundary conditions far exceeds the lateral curvature term on the finger boundary. 
The reason for this unexpected validity of the MS theory is, first, that both Ca and 
Q are small so that the deviation from the Saffman-Taylor finger solutions is actually small. 
Indeed, the role of terms such as lateral curvature and transverse curvature is not so much 
to change the Saffman-Taylor shapes as to determine the finger width which is arbitrary to 
zeroth order. Second, the finger width is being determined by transcendentally small terms 
in Q in the physical domain which can only be determined by analytic continuation of the 
equations to the neighborhood of some point in the unphysical plane that is the source of 
the transcendentally small correction. The relative size of lateral, transverse curvature and 
thin film leakage terms in an ’inner’ region near this point is rather different from what 
they are in the physical domain; yet it is this relative size which determines the finger 
width. It turns out that the lateral curvature is far bigger in the ’inner’ region than the 
thin film leakage or the transverse curvature term when Ca << Q 7/ ^ 3 << 1 and that 
explains the unexpected persistence of the solutions based on the MS theory. However, 
in this paper, for this range of parameters, we do not address the question of existence of 
other kinds of solutions not found in the MS theory. 

Our analysis also shows the importance of doing a nonlinear analysis in the region 
of nonuniformity of the original perturbation expansion. We point out that in the case 
of the McLean-Saffman boundary condition, nonlinear analysis is only necessary on an 
equation containing one parameter . If the nonlinear equation is replaced by 

a linear equation, the linear equation still contains this parameter. Only the numerical 
value of • obtained by matching to the outer solution will be affected by the ad 

hoc linearization. Thus, only the proportionality constant between A — ^ and C? 2 ^ is 
affected if one resorts to a simpler linearized analysis. With the thinfilm effects, there 
are, in general, nested inner regions surrounding the point of non-uniformity of the outer 
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perturbation expansion and the form of the nonlinear equations together with the type of 
parameters in each such region is determined by the nonlinearity of the equation in the 
next outer region. If the equations were linearized at the outset, one is unable to identify 
the correct scaling law between the three parameters A , Q and Ca , let alone determine 
the scaling constants. 

2. Mathematical Formulation and leading order solution: 

We consider the two-dimensional averaged flow in a Hele-Shaw cell of width 2 a where 
a finger of zero viscosity, moving with velocity U , displaces a viscous fluid of viscosity n . 
The 2-D averaged velocity field in the region not occupied by the finger is the gradient of 
a harmonic potential function (f> . The finger boundary is assumed to be asymptotically 
parallel to the walls. Here, in this paper, we will take each of a and U to be unity without 
any loss of generality since this is equivalent to non-dimensionalizing all our variables using 
these two parameters. We will now proceed with the understanding that all the variables 
have been non-dimensionalized. The flow domain in the frame of the steady finger, the 
z = x + iy plane, is shown in Fig. 1. We introduce the complex velocity potential, 
W = <f> + i ip , which will be an analytic function of the complex variable z in our flow 
domain. In the frame of the steadily moving finger, the boundary conditions on the finger 
that incorporate the the transverse thin film effect (see Reinelt (1987b) for details) are: 


, , 4 g 

q> + x = 

7T Z 


( 2 . 1 ) 


= 


m°(Ca n x ) + m 1 (Ca n x )\J ~ £/ 1/2 Ca 1/2 i 


( 2 . 2 ) 


where Q = , Ca = ^jr (capillary number), n x = cos 6 is the component of 

the outward unit normal to the finger (see Fig. 1) along the x -axis, u is some parameter 
in the interval (0, 7 r) (specific choice given later) parametrizing the finger boundary and 
the subscript with respect to u denotes derivative with respect to u . R is the radius of 
curvature of the interface in the lateral ( x — y ) plane. The functions , m° and 
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m 1 have have been calculated by Reinelt (1987a) numerically over a range of arguments. 
The appearance of n x in the argument of these functions is because only the normal 
component of the finger velocity is relevant in the thin film effect. In the limit of small 
Ca , following Bretherton, Reinelt found that 


K°(Ca) = constant — 3.878 Ca 2 / 3 + ... 

(2.3) 

K 1 (Ca) = — 7r / 4 + 4.153 Ca 2/3 + ... 

(2.4) 

m°(Ca) = 1.3375 Ca 2/3 + ... 

(2.5) 

m 1 {Ca) = -1.3375 ^ Ca 2/3 

(2.6) 


The values of the constant coefficients in (2.3) and (2.4) are different in their third signifi- 
cant figures from those originally quoted in Reinelt (1987b). Reinelt (private communica- 
tion) found these revised values after a more accurate calculation than previously reported. 
Note that the constant term in (2.3) does not affect the finger shape or the fluid velocity 
field since it can be always absorbed as part of x on the left hand side of (2.1) by suitably 
choosing the origin of x . Thus from this point onwards, we will neglect the constant term 
on the right hand side of (2.3). The MS theory corresponds to neglecting all other terms 
involving powers of Ca , while replacing by — 1 . A rational approximation when 
Ca << 1 and Ca << Q 3 would be to replace /c 1 by — ^ while neglecting all terms 
in (2.1) and (2.2). However, this only modifies the sole control parameter Q of the MS 
theory by a factor of \ . Also the range of Ca and Q investigated in the experiments 
to date, it is clear that Ca is not small enough to justify this approximation. 

The mathematical formulation here is a generalization of the previous formulation 
(Tanveer, 1987b) developed for the MS boundary conditions. In the previous work with 
the MS boundary conditions, no assumption was made on the symmetry of the finger about 
the channel centerline. As part of the conclusion, it was found that only symmetric fingers 
are possible for non zero Q . In this case, for the sake of simplicity, we will assume apriori 
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that the finger is symmetric about the channel centerline. The existence of nonsymmetric 
finger with the SPHR boundary conditions is still an open question. 

Consider the conformal map of the flow domain in the z -plane into the interior of 
the unit semi-circle in the C -plane (Fig. 2) such that z = oo, z = — oo + i and 
z = — oo — i are mapped to C = 0 , +1 and -1, respectively. The finger boundary 
is mapped to the arc of the semi-circle in the C -plane. Introduce analytic functions 
/(C) and g(() in the unit ( -semicircle so that 

* = -- In C + - (1 - A) In (C 2 - 1) - i (1 - 2A) + - /(C) (2.7) 

7T 7T 7T 

IF(C) = —2 ^— In (C 2 - 1) + 2(1 ~ A + a) InC + i(l-A-a) + - g(C) (2.8) 

7T 7T 7 r 

where 

a = — f du m° + m 1 J ^ 1/2 Ca 1 ^ 2 y v dv (2.9) 

L T 7T Z itj 

where C = e' v on the semi-circular boundary. Note a accounts for the total flow across 
the interface edge. From conservation of fluid flow, it is easy to show that the fluid velocity 
V far ahead of the finger in the laboratory frame is given by V = (A — a) U . Then the 
condition that the real diameter of the unit semi-circle in the C plane corresponds to the 
walls of the cell where Im z = ±1 implies 

Irn / = 0 (2-10) 

The condition that on the upper and lower cell walls Im W = =p(l — A + a) implies 

Im g = 0 (2.11) 

The boundary conditions (2.1) and (2.2) can then be written as the following boundary 
conditions on / and g on the semi-circular boundary Q = e xv , where v is in the 
interval (0, 7r) : 

Re (/ + g) = - 3~ 1/2 Q 1 ' 2 Ca~ Xl2 k° - — ^ (2.12) 

7T R 
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Re(Cg') = -m° Re C (/' + fc) - m 1 )/ £ 1/2 Ca 1/2 Re C (/' + *0 - a ( 2 - 13 ) 

where primes denote derivatives with respect to the argument and 


MO 


(i-p 2 C 2 ) 

C(C 2 -1) 


(2.14) 


where 

p 2 = 2A — 1 


is in the open interval (-1, 1) for A in the open interval (0, 1). Note that a can be 


expressed as: 


a 



m 


+ m 



Re e iv [f'{e iv ) + h (c* 0] (2-15) 


Also note that owing to Schwartz’s reflection principle each of (2.12) and (2.13) hold in 
the lower half semi-circular boundary as well. 

One can establish that on £ = e iV , u in the interval (0, 7 r) , 


n x 


1 

R 


l r + M 


Re [C (/' + h) ] 


TT 


2\(f' + h)\ 


Re 


l + c 


r + h r 

f' + h ■ 


(2.16) 

(2.17) 


There is an alternate way of writing each of (2.16) and (2.17) which defines analytic 
continuation of n x (0 and ^ (C) °ff the arc of the upper half unit £ semi-circle. 


where 


MO 


1 hiQ ~ *i(0 

2 (/l(0) 1/2 (^2(C)) 1/2 


R® 


7ri(C 2 — 1) l 1 + 


+ 


1 

5 




(ii(0) 1/2 (WO) 


1/2 


h( 0 = C 2 -p 2 + i/'( l/O 
h( 0 = 1-P 2 C 2 + C(C 2 -i)/'(0 


(2.18) 


(2.19) 


( 2 . 20 ) 

( 2 . 21 ) 
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For Ca = 0(1) and Q << 1 , the leading order approximation can be found by 
setting Q = 0 in (2.12) and (2.13). Thus, to the leading order, g = — /q , / = /o , 
where fo is determined by the boundary condition 

Im f 0 = 0 (2.22) 

on the diameter of the unit semi-circle and on the arc of the unit circle, £ = e iv 

R e C fo = m {Ca n x(C) ) Ra [C (fo + ^)] + (2.23) 

where is n z with the substitution f = fo in (2.16) and a 0 is the constant a as 
determined by (2.15) when m 1 is neglected and Can j? is substituted as the argument of 
m° . The solution to (2.23) was found numerically for arbitrary p 2 in the open interval 
( 1 * 1) > A in (0,1) , by expanding /o in a power series 

fo = a l C 2 + C 4 + •• (2.24) 

truncating the series to N terms and satisfying (2.23) at N uniformly spaced out points 
on the arc of the unit ( quarter circle. The condition (2.10) is automatically satisfied by 
restricting to real a® . The resulting system of nonlinear algebraic equations was solved 
using Newton iteration and consistency of results checked by doubling N . 

Higher order formal corrections to the perturbation expansion in powers of Q ^ 2 can 
be found by expanding / and g as: 


/ = fo + G 1 ' 2 h + .. 

(2.25) 

4 

9 = 90 + G 1/2 91 + •• 

(2.26) 



where fo and go are as determined above. f\ and g\ are determined by satisfying 
equations: 

Im f i = 0 (2.27) 

Im g\ = 0 (2.28) 
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on the real diameter and on the arc of the unit semi-circle in the £ -plane, 


Re (/i + £i) = -3 ~ 112 Ca~ l/2 K°{n° x Ca) (2.29) 

Re £ g'i + m° (Co n°) Re (f[ + m°'(Ca n°) Ca Re £(/q + h) n\ 

= -m\Ca n°) \j H Ca 1 ' 2 T- fie C(/o +*)-“! ( 2 - 30 ) 


where 


n, = — 


1/6 + *1 


fie C/i + 


fie C(/o + ft) 

1/6 + M 


Re 


(_iU 

'/i + *' 


For each A and Co , by truncating power series representations in powers of £ 2 for each 
of /1 and <71 to a finite number of terms and satisfying (2.29) and (2.30) at a discrete 


set of uniformly spaced out points on the unit £ quarter circle, we solved the resulting 


system of linear algebraic equations. The resulting a n and b n were found to decay with 

n for large n and there was consistency on doubling the collocation points and unknowns. 
This suggests that a regular perturbation expansion fails to select A in contradiction to 
direct numerical results of Reinelt (1987b). A similar contradiction occurred earlier in the 
context of McLean-Saffman boundary conditions (McLean & Saffman, 1981) which was 


later resolved (Combescot et al (1986), Shraiman (1986) and Hong & Langer (1986)) by 
extracting transcendentally small terms in Q . 

When Ca « 1 , the absence of selection from terms of a regular perturbation 

expansion can be more clearly demonstrated. Since the right hand side of each of (2.12) 
and (2.13) are small for small Ca and Q as is clear from (2.3)-(2.6), to the leading 
order, / = 0 and g = 0 , i.e. we get the Saffman Taylor finger solutions in which 
A and hence p 2 is arbitrary. Now, if we were to try to improve this approximation by 
including the next order correction, we would substitute / = 0 and g = 0 into 

each term on the right hand side of (2.12) and (2.13) and determine the first non trivial 
g and / term on the left hand side. One can systematically proceed to determine each of 
/ and g in a perturbation expansion involving powers of G 1/2 and Ca 1/6 , the precise 
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form depending on the relative ordering of the two small parameters. This can clearly 
be done again without any constraint on p 2 , since solving at each stage corresponds to 
calculating two harmonic functions Re (/„ + g n ) and Re g n with specified Dirichlet 
and Neumann data respectively on the unit circular boundary where the Neumann data 
consistency condition is easily seen to be satisfied. The solution g n is only determined 
to an arbitrary additive constant that does not affect fluid flow or the finger shape. There 
are no sources of nonuniformity of this outer perturbation expansion except possibly near 
C = ±1 > i-e. the tail of the finger. For the MS boundary conditions, McLean & Saffman 
(1981) have shown that a secondary expansion near the tail is possible and that this 
matches with the outer expansion without any constraint on A . We expect that such a 
matching is possible in this problem as well. We will assume this is indeed the case. Thus, 
given this assumption, it is clear that A remains undetermined by a regular perturbation 
series. 

Thus, to calculate A analytically for any Ca , one must calculate terms for / and 
g that are transcendentally small in Q . These are terms beyond all orders of the regular 
perturbation expansion in powers of Q and are subdominant except when every term of 
the regular perturbation expansion is zero. Kruskal & Segur (1986), in their pioneering 
work in calculating terms beyond all orders for a model third order nonlinear ODE, give 
a careful account of the issues involved in calculating terms beyond all orders. 

As with the earlier case of MS boundary conditions, it will turn out that there is 
some open interval on the imaginary C axis containing i where the condition of a 
finger boundary with a smooth tip, Im f = 0 , is satisfied by every term of the 

regular perturbation series of / in powers of Q^^ 2 for any A ; yet the leading order 
transcendentally small term in Q violates this condition except when A is restricted to 
a discrete set of values. 

We note that to find transcendentally small terms in the physical domain, we must 
find the source of non-uniformity of the regular perturbation expansion in the unphysi- 
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cal domain. Such sources of non-uniformity of the perturbation expansion contribute to 
transcendentally small terms in the physical domain. One needs to rescale dependent and 
independent variables in the immediate neighborhood of the source of non-uniformity and 
construct an inner perturbation expansion. An inner-outer matching is then carried out 
and the terms that are beyond all orders in the matching procedure are transcendentally 
small in the outer region that includes the physical domain. In this paper, only the leading 
order inner and outer matching is carried out since only the leading order transcendentally 
small terms will be found. As a first step, it is necessary to locate the source of nonuni- 
formity of the regular perturbation expansion and construct an inner expansion near the 
source of non- uniformity. For that purpose, it is necessary to continue analytically each of 


(2.12) and (2.13) outside the unit C circle across the arc of the upper half semi-circle. 

3 . Analytical continuation of the equations to |£| > 1 

First, we note that from Poisson’s integral formula relating a harmonic function and 
its conjugate inside the unit circle to its value on the boundary, one finds that for any 
analytic function -F(C) in the unit semi-circle with vanishing imaginary part on the real 
diameter, 


F(0 


2n i 


f dC 

(C + C 

1 + CC'l 

1 ft 

c s 

k'-c + 

l - CC' ■> 


(3.1) 


where the contour of integration C in the (' plane is along the the semi-circular arc 
from +1 to -1. For convenience, we define the operation on Re F on the right hand side 
of (3.1) as l(Re F) . If Re -F(C') = G(C') for £' on the unit semi-circular arc, where 
G(C) is an analytic function defined off the semi-circular arc as well, then by deforming 
the contour in the (' plane, one finds that the analytic continuation of F( C) across the 
semi-circular arc for |£| > 1 must be 


F(C) = 1(G) + 2 G(Q (3.2) 

Note that on the arc of the unit circle, £'* = £ _1 and Re F = + -F(C -1 )] f° r 

any function F which has vanishing imaginary part on the real diameter of the unit 
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semi-circle. Using this and the continuation property (3.2), one finds that the analytic 
continuation of (2.12) and (2.13) to |C| > 1 across the arc of the upper half semi-circle 
must be 


/ + 9 — MC) ~ 


4 Q o 


3 Ca 


k (Ca n x ( 0) + /12(C) 


-7 4(0 K 1 (Can x ( 0) (3.3) 

7T it 


( s' = MC) - m ° ( Ca n «(0) + MO 

- 5<0 h(0 (2 -_ I 110 - l (Ca MO) - « (3-4) 

where we assume each of the functions m° , m 1 , k° and k 1 are locally analytic 
functions of their argument and 


MO = -v ^ I ( K°) (3.5) 

MO = - v 1 ( 4) (3.6) 

7 r K 

/13(c) = ( m ° ^ 2 _ j 1 ) ( 3 - 7 ) 

MC) = -yj I (rn 1 i f ^ 2 _ j 1 ) (3.8) 

Note that each of the functions h\ , fi2 , /13 and /14 involve integration of values of 
functions / and g on the arc of the unit ( '/ circle and therefore analytic everywhere in 

ICI > 1 • 

Equations (3.3) and (3.4) are a pair of nonlinear integro differential equations for 
/ and g that appears too formidable to be of any practical utility. However, it is 
amenable to simplifications when Q is small. When |C| > 1 , |C _1 | < 1 and so terms 
such as /'(1/C) , /"(1/C) appearing in j and n x in (2.18) and (2.19) can be replaced 
by the regular perturbation (2.25) even when the deviations of /(C) from the regular 
perturbation series (2.25) is not small as is true near the point of nonuniformity of the 
regular perturbation series. To the leading order, it is appropriate to replace f(l/() by 
/o(l/C) and /"(1/C) by fo (1/C) • The same substitution is appropriate for / appearing 
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in the integrand in each of h\ through h 4 . Then the pair of equations (3.3) and (3.4) 
become a pair of nonlinear second order differential equation for / and g . We are 
interested in solutions to these equations that match with (2.25) and (2.26) as one moves 
away from the points of nonuniformity of the regular perturbation expansion towards the 
arc of the upper half unit semi-circle. 

4. Nature of transcendental correction 

We now seek the form of the leading order transcendentally small correction terms 
in Q to the regular expansion (2.25) and (2.26). As has been discussed before (see, for 
example, Kruskal & Segur (1986), Combescot et al (1986), Tanveer (1987b), Combescot 
et al (1988), Combescot & Dombre (1988) and Tanveer (1989)) linearizing the nonlinear 
equations about the leading order behavior given in (2.25) and (2.26) and looking for WKB 
type of solutions to the homogeneous part of the linearized equation. 

In our case, the resulting homogeneous equations arising by linearizing (3.3) and (3.4) 
about / = / 0 + fi are: 


fa + 9 h - Q 1/2 f' H ~ j- [ I'm + Q 2 J'h ] 

C q'h = C Qz /h + C G 1/2 [ J'h + Q 2 I'h ] 

where 


Q 1 


^2„ + lu 
4« 1 C(C 2 -1) 



k 0 ' ( Ca n°) + 


4^72 

7T Rq 


Ca k 1 ' (Ca n^) 


(4.1) 

(4.2) 


(4.3) 


L = 


-HU 1 ' 2 *2, 3/2 


k 1 (Ca n£(C)) C 2 (C 2 - 1) 

0 = 3 C 2 -1 _ £ !k , I _ I % (VO 
V2 c (C 2 - 1) 2 l 2o c 2 h 0 


(4.4) 

(4.5) 


Q 3 = - m°(Ca n° t (()) 


1 12 Q C a 1 1 . n , 

- V -72— * 0 "* ( Cu »- (0) 


- Ca l ^3/2^1/V 02. -' 1 .) {™°' (CanJ) + \j 12 y n 2- mV {Ca "■ <0)_ 

2o lo 


(4.6) 
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Q 4 = - V 


J 12 Ca l 2o - Zl 0 m\Can ° (0) 


i?0 


16 C(C 2 - 1) K 1 (Ca n«(0) 

a f 1 I 1^0 (0 I 1 ^3 0 (VC)l 


(4.7) 


(4.8) 


2 (U0) 1/2 (MC)) 1/2 

where the zero subscript on each term here refers to the evaluation of that term with the 
the subsitution / = /o + (7^ 2 /l ■ Eliminating g H between (4.1) and (4.2), one obtains 
the following second order linear homogeneous equation for f' H 

Q 4 


in + IQa + jaS + $ 71 ~t) % 


+ /» ( 

In the limit of Q 


n' L> n 4 . Q i 

Q 2 ~ 1 Qi+gin 


QH 2 ' & ] 

Qi L 1 L . <?4 n 

2 + TrU + ^3) + 7TT79 V2 


L Ql> 2 ' ' Q 1/2 

0 , one can seek WKB solution of the type 


)- 


0 (4.9) 


f' _ W 0 + w, +o(c 1/s ) 

Jff ~ e 


(4.10) 


to find that 


W/, = 


Q4 + Q 1 


«?4 + Ql)" 


- L (l + Qs) 


(4.11) 


W{ = 


■WS - wft (g 2 - ^ ) - Q 1 + Qi t ~ Q4Q2 

2 IEq + Q 4 -(- Qi 


(4.12) 


Let us define g\ and g 2 as the asymptotic solution in (4.10) corresponding to the two 
different choices of the sign of the squareroot in (4.11). Thus 


and 


9i 


p G~ l/3 W+ + W* 


(4.13) 


9 2 


,-i/j 


= e” 


w ' 


+ w r 


(4.14) 


where the plus and minus superscript on Wo and W\ correspond to the choice of plus 
or minus sign in (4.11). Further, to the order of approxmation that the WKB solutions 
(4.13) and (4.14) are valid, it appropriate to neglect £/ 1/2 in the zero subscripted terms 
in Qi , Q 2 , Q 3 , Q4 and L appearing in the right hand side of (4.12) but not (4.11). 
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The condition Re = 0 and Re Wq = 0 (caUed the Stokes lines) divide 

the upper half C plane into different sectors. Including the leading order transcendental 
correction in the description of the asymptotic behavior of / in some sector of the complex 
£ plane adjoining the arc of the upper half unit semi-circle |C| — 1 , Ivn C > 0 , we 
can write 

/ = /o + G 1/2 fl + ••• + Cl gi + C2<72 ( 4 - 15 ) 

where the transcendental term C\ gi + C 2 92 is assumed small. It is generally sub- 
dominant to every other term of the perturbation solution denoted symbolically by ... in 
the above equation. If Re >0 gx is transcendentaUy large. If Re Wq > 0 , 

g 2 is transcendentaUy large. In each case, the corresponding coefficient C\ or C 2 has 
to be zero since each of Re and Re Wq increase without bounds as ( = ±1 or 
C = 0 is approached as can be seen from the integration of (4.11) once expressions for 
Ql , <54 and L are substituted. As we cross any Stokes line, the coefficients C\ or 
C 2 in (4.15) can change. The actual determination of the coefficients C\ and C 2 is 
possible only by matching an inner perturbation expansion carried out in the neighbor- 
hood of points where the outer perturbation expansion (2.25) and (2.26) to (3.3) and (3.4) 
as well as the leading order WKB solutions (4.13) and (4.14) to (4.9) are invalid. This 
is rather an involved procedure for arbitrary Ca since the exact solution for /o is not 
known and even the source of nonuniformity of the outer perturbation expansion has to be 
determined by a nontrivial numerical procedure. Further, no analytical formulae for the 
functions ro° , m 1 , k° and k 1 exist for arbitrary values of their argument. In section 
7 , we will show more precisely what pieces of information would complete the selection 
theory for arbitrary Ca and give indications on how they might be obtained. 

From now on till section 7, we will limit ourselves to small Ca and assume the 
formulae (2.3)-(2.6) to be valid. In this case, as noted earlier, the Saffman- Taylor exact 
finger solution / = 0 and 5 = 0 provides a convenient starting point. Higher order 
approximations involve powers of Ca l/6 and G 1/2 , the precise form depending on the 
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relative ordering of the two small parameters Ca and Q . We will symbolically denote 
the leading order correction with algebraic dependence on Ca and Q by fi . Thus, away 

from the points of nonuniformity of the regular perturbation expansion, which in this case 
are at ( = ±- , 

/ ~ fl (4.16) 


To this, one can add terms with higher order algebraic dependence on Ca and Q . The 
behavior of the analytic continuation of f\ in the region |£| > 1 can be obtained from 
(3.3) and (3.4) by setting /' and f" to zero in every term involving these on the right 
hand side of (3.3) and (3.4) after the Bretherton approximations (2.3)-(2.6) are invoked. 
One can see that /i is singular at £ = ±1 /p which are then the sources of nonuniformity 
of the regular perturbation expansion. We also r f ice that f\ is real on some open interval 
on the imaginary axis that includes £ = * • The same can be shown to be true to every 
order of the regular perturbation expansion. However, transcendentally small corrections 
to (4.16), as we shall see in a moment, do not generally satisfy this condition implying 
that the finger tip is not smooth for arbitrary A . 


Terms that are transcendentally small in Q , which to the leading order are constant 
multiples of g\ or g 2 as determined by (4.13) and (4.14), can be added to (4.16). When 
|1 — p 2 £ 2 | >> Ca 2/7 -j , the expression for the asymptotic WKB solutions to (4.9) 

simplify even further because the capillary number Ca has been assumed to be small and 


so each of Qi , Q 2 , Qz and Q 4 can be neglected. In that region, to the same order 
that (4.13) and (4.14) are asymptotic solutions, g\ and g' 2 can replaced by g{ and 


g 2 where g\ and g 2 are defined as 


91,2 


a -P 2 ef ,s 

(C 2 - p 2 ) } 


± 


e 


i 3 / 3 

fli/j 



d( 


Ci- y 3 C 3 ) 3/ * 




U a - V 7 ) 1 '* 

- o 


(4.17) 


Note that g\ and <72 <M:e transcendentally small or large depending on the sign of Re P , 
where 


P = 


,3/2 




dc 


(1 -P 2 C 2 ) 3/4 (C 2 - P 2 ) 1/4 
C (C 2 - 1) 


(4.18) 
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The Stokes lines here are the same as for the M clean- SafFm an boundary conditions as dis- 
cussed earlier (Tanveer, 1987b). We have different cases depending on p 2 < 0 or p 2 > 0 . 
The Stokes lines in each case are shown in Figs. 3 and 4. Note that p 2 < 0 corresponds 
to A < \ and p 2 > 0 corresponds to A > ^ . 

5. Restriction to A < 5 

It is appropriate in this case to define q , a positive real parameter, so that p — iq . 
In that case the point of nonuniformity of the regular perturbation expansion in the upper 
half plane £ = jr = f as shown in Fig. 3. 

Away from the immediate neighborhood of £ = ^ > the desired asymptotic behavior 
for / must be 

/ = fi + •• + Cl 91 (5-1-) 

in sector I of Fig. 3; in sector II, one must have 

/ = fi + ... + C2 92 (5-2) 

and in sector III 

f — f l + .... + C3 52 ( 5 - 3 ) 

because g\ is transcendentally small in sector I, and 52 transcendentally small in sectors 
II and III. No transcendentally large terms are allowed in this three sectors as each of 
the sectors I, II and III extend all the way to the physical region inside the unit upper 
half £ semi-circle. If the finger boundary is smooth, the solution / to the analytically 
continued equations (3.3) and (3.4) that match with the regular perturbation expansion 
as the physical region is approached must be unique and thus (5.1), (5.2) and (5.3) are the 
different behavior of the same global analytic function continued across the Stokes lines. 
Further, we note that 51 and 52 are both real on the imaginary axis in the interval 
(j 5, -) and therefore, the assumed finger symmetry about the chnnnel centerline implies 
C2 = C3 . Thus, to the leading order, it is enough to require that each of (5.1) and (5.2) 
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hold in sectors I and II and in addition require 


Im C i = 0 (5.4) 

Generally, for any A , it is possible to find asymptotic solutions to (3.3) and (3.4) 
satisfying (5.1) and (5.2). However, (5.4) is violated. This implies that the tip is not 
smooth. Further, from the behavior Im f ~ Im C\ g\ on the imaginary axis in the 
vicinity of £ = i , one is able to deduce that the singularity of / at £ = i is a branch 
cut corresponding to a cusp at the finger tip. At this point, it is appropriate to point out 
that this generalized solution with cusp is an asymptotic solutions to (3.3) and (3.4) for 
|C I > 1 in the 1st quadrant that matches with the perturbation series (2.25) and (2.26) 
as the physical region |£| = 1 is approached from the exterior. It is not an asymptotic 
solution to the original equations (2.12) and (2.13) on the entire arc of the entire semi- 
circle as a cusp would imply infinite curvature at the tip which is clearly not possible if 
(2.12) were to be satisfied at the tip as well. A solution to (3.3) and (3.4) is also a solution 
to (2.12) and (2.13) only when we require a unique analytic continuation of / across 
the semi circular arc satisfying conditions (5.1), (5.2) and (5.3); which for a symmetric 
finger is equivalent to satisfying just (5.1) and (5.2) along with the requirement that on 
the imaginary ( axis there is some interval containing ( = i where Im f = 0 . To the 
leading order, this latter condition is equivalent to (5.4). To find the actual expression for 
C i , it is necessary to find and solve inner equations near £ = i where both the outer 
perturbation expansion and the WKB solutions g\ and <72 are invalid. In the following, 
A together with Ca and Q will be treated as parameters and solutions found that satisfy 
(5.1) and (5.2). Once such solutions are found, the condition (5.4) will be imposed to see 
if they can be satisfied for some constraint relation between the three parameters. 

The details of the inner region for A < ^ depend on the assumptions made on Ca , 
Q and A . Here w r e will consider three different cases: 

(a) (i_ 2 xyQ << 1 or order unit y> 

(b) (1 — 2A) 2 C/A -4 << Ca «G 1 ' 2 ( l-2A)A- 2 and 
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>> 1 and in 


(c) Ca Q ~ ^ 2 (1 — 2A) -1 A 2 >> 1 or order unity. 

It will be assumed that in case (a), r >> 1 , in case (b) r /?{ 3/2 
case (c), r >> 1 , where r , f3\ and are as in (5.8), (5.11) and (5.12) 

below. However, as shall be seen later, a sufficient but perhaps not necessary condition 
for the above assumptions to be valid is that Ca = 0(Q 7 ^ 3 ) or larger in case (a) and 
Ca = 0(G ) in case (b). In case (c), the necessary assumption is equivalent to requiring 
that Ca -7 / 4 A" 7 ' 2 (1 — 2A) 7 ^ 2 >> 1 . The solutions found in this case are found to 

satisfy this condition aposteriori provided Ca >> Q . However, the possibility is open 
for new kind of solutions not described here with A < 1/2 but approaching ^ in the 
limit Q — > 0 provided Ca = 0( Q) or smaller. 

5a. Case o f (i_ a -2\yg << 1 or order unity. 

In this case, it is appropriate to introduce introduce inner variables 


/ = 


2T 

l + q : 


F(t)r 


- 4/7 


(5.5) 


9 = 


2<r 


o>(0 r 


1 + q 2 
1 + iqC, = r- 2/7 £ 


- 4/7 


where 


r = 


2 3 / V(1 - 9 2 ) 1/2 4 


(5.6) 

(5.7) 

(5.8) 


(1 + q 2 ) 3,2 Gn 

We will assume that (1 — 2A) 3 ^ 2 A 1 / 2 >> G so that r is indeed a large parameter. 
This assumption can become invalid in two ways: first, if A is so small that A ~ C/ 2 , 
and secondly when A ~ ^ such that 2 , ^ = 0(1) . For the first case, since 
A >> e , i.e. A >> G^ /2 Ca}! 2 for the SPHR boundary conditions to be valid, it is 
clear that the assumption can possibly be violated only when G ^ 2 Ca 1/2 << G 2 , i.e. 
when Ca << Q 2 • However, a glance at the original equations (1) and (2) would seem 
to suggest that under these conditions this film effects can be neglected and that the MS 
boundary conditions are valid. We will assume this is indeed the case. Earlier work on the 
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MS boundary conditions rule out any possibility of A < 1/2 . For the second case when 
A ~ ^ , the assumptions made in this case imply that (1 — 2A) = 0(Ca ^ 2 Q~ or 
larger. Thus the assumption of r large can be violated only if Ca << Q 7 ^ 3 . Thus if 
Ca = 0( g 7/3 ) or larger, the assumption of large r follows. However, this may not be 
a necessary assumption. When r is not large because q is small, one can use the scaling 
in section 6 with p replaced by q ; however the analysis for such a case is not presented 


here. 


From the transformations (5.5) and (5.6), to the leading order in r 2 ^ , we find that 
(3.3) and (3.4) reduces to 


and 


where 


F + u> = -0i 


u/ = 02 


+ 


(F» - 1) 
(£- F ') 1/3 ’ (£- F ') 3/2 

(F" - 1) 


+ 03- 


01 = 


3.878 ^ Ca^ 1/6 (1 + tf 2 ) 4/3 (1 - g 2 ) 1/3 


2 \/3 V Q 


qS/3 r -2/3 g- 2/3 


a _ i ooyr Cfl2/3 ( X ~ , q^/21 

02 — 1.3375 2 i A S/3 {r y) 


Q*l 


4 ^ 8/3 


and 


03 = 1-3375 


3 7r Ca 


7/6 


(1 - O 


4 . 5/6 (1 + q 2 ) , ,^17/21 


2 32 ^13/ 42 g 

Eliminating w between (5.9) and (5.10), we obtain 


11/3 


(r ^) J 


(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 


F' + 


02 


+ /?3 


(F" - 1) _ 0i (1 -F») 


+ 


3 (F" - 1)' 


({ _ pj 1 ' 3 ' (( - F’J 11/G 3({ - F') 4 ' 3 ({-F') 3 ' 2 2 (« - F') 5 ' 2 

(5.14) 


Note that in order to obtain (5.9) and (5.10) from (3.3) and (3.4), we have to assume 
Ca << £/ 1/7 A - 4/7 (1 — 2A) 2/? or otherwise the product Ca n*(C) f° r £ = G(l) is 
not small making the Bretherton approximation (eqns 2. 3-2. 6) invalid. 

Further, in this case 0\ , 02 and 03 are far smaller than unity and one can neglect 
to the leading order all terms in (5.14) involving 0\ , 02 and 03 . Going back to (5.9) 
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and (5.10), this means that to the leading order oj = 0 without any loss of generality 
and 


F = 


(F" - 1) 


(5.15) 


(£ _ F'f 2 

which is the same equation as for the MS case when thin films and the variation of trans- 
verse curvature are totally ignored. The solution F to (5.15) matches with the outer 
solution in sectors I and II provided we require that the F goes to zero as £ — ► oo for 
Arg £ in the interval [0, 67t/7) . This requirement gives a unique solution to (5.15) and 
thus matching with (5.1) and (5.2) for |l + *'?C| << 1 , one finds C\ and C 2 . The details 
of this case have been discussed in the context of the analytic theory of velocity selection 


of a bubble using the MS boundary conditions (Tanveer, 1989a) and also in the context 
of dendritic crystal growth (Tanveer, 1989b), where the same equation (5.15) appears. 
Im C\ can only be zero if the solution to (5.15) satisfying the above stated conditions is 
real on the real axis and this is found not to be the case. Thus there are no solutions in 
this case. 

When Ca = 0(Q (1 — 2A) 2 A -4 ) , f3i = 0(1) , while other parameters and 
ft are small. Thus to the leading order, one can ignore the j3 2 and ft terms in (5.14). 
This means that in this situation, we can ignore the thin film effect on the kinematic 
equation but include it in the dynamic conditions. Then going back to (5.9), the relevant 
equation in this case to the leading order is 


F 


= -ft 


a ~ F') v 3 


+ 


( [F " - 1 ) 
(£ - F'f /2 


(5.16) 


Numerical solutions of this equation for /?i in the range (0, 10) were sought and a 
unique solution to (5.16) was found that tends to 0 as ( -► oo for Arg £ in the interval 
[0 67T / 7) as necessary to match with (5.1) and (5.2). However, such a solution was not real 
on the real £ axis as would have to be the case to satisfy (5.4). The numerical method 
is a trivial modification of the earlier procedure (Tanveer, 1989b) for finding appropriate 
solutions to (5.15). Indeed with increase in , the residual of the numerically imposed 
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condition for a smooth tip increased with increasing pi suggesting that an increase in 
Pi causes larger mismatch in the tip slope between the two sides of the cusp. 

5b. Case of (l-2A) 2 £A~ 4 << Ca « G 1/2 ( l-2\)\~ 2 

In this case, Pi is large. pz is small; but P 2 can be large or small compared to 
unity. It is appropriate to introduce rescaled variables 


6 = ft’ 3 ' 7 r 2 ' 7 (1 + ; q C) (5.17) 

ft (Cl) = -ft- 6 ' 7 r 4 ' 7 /(C) (5.18) 

u-i(fl) = -ft' 6 ' 7 ^p-r 4 ' 7 s (() (5.19) 


where r and Pi are as defined earlier in (5.8) and (5.11). 

We will assume that r py 3 ^ 2 >> 1 , implying Ca~ *^ 4 G~ 3 ^ A" (1 — 2A) >> 1. 

The only way this assumption can be violated is if A is very close to ^ • However, 
in this case, when A ~ ^ , (1 — 2A) >> Ca G~ 1/2 . Thus a sufficient con- 
dition for the validity of the assumption is that Ca 0( G) or larger. Note that the 
Bretherton approximation is valid, i.e. Ca n x << 1 for = 0(1) , since 

Ca n, = 0(Ca 27/28 Q~ 3/28 A 3/7 (1 - 2A)~ 3/14 ) << 1 . This is true because 


(1 — 2A) » and Ca 3 / 4 « 1. 

Then to the leading order in p 3/7 r~ 2/7 , for = 0(1) , equations (3.3) and (3.4) 
become 


u>l + Fi = 


, (h 

Wl = T477 


1 , a- 3/2 gj 

(Ci - ft) 1 ' 3 1 (ci - F o 3 ' 2 

ft (ft' ~ 1) 


W ' 7 (Ci - ft ) 1 ' 3 + ft 17 ' 14 (Ci - n ) 11 ' 6 


(5.20) 


(5.21) 


-3/2 


From the definition of Pi , P 2 and Pz , it is clear that each of the quantities Pi , 

-Att and are much smaller than one. Thus to the leading order, the appropiate 

f>\ * ••• - ' ^ 

equation is „ 


ft. = 


(Cl - Fi) 1 ' 3 


(5.22) 
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For large £1 , the solution must match with / ~ f\ . On examination of (5.22) the 
behavior of -Fi(£i) needed to assure this matching is: This implies that for large £1 , the 
asymptotic behavior of F\ 0 must be 


Fi b 


d ' 3 


(5.23) 


On linearizing (5.22) about this behavior, we obtain the homogenous equation 


F' h + 3 £l 4/3 

The solution to this is of the form 


9^ /3 


F B = 0 


(5.24) 


F Hl = ft 4/3 e~H 


• t v» 


(5.25) 


This is transcendentally small for Arg in (—3 7r/14, 37 t/ 14) and transcendentally 
large for Arg ( i in (3 tt/14, 97t / 14) . A second possible transcendental behavior for large 
is found by retaining the highest derivative term in (5.20). Linearizing about (5.23), 
and examining the homogeneous equation it is easy to see that a possible transcendental 
behavior for large £ is 

F H3 = fV' f «: /# [l + o(l)] (5.26) 

and this is transcendentally large for Arg in the interval (— 37r/7, 3n/7) and tran- 
scendentally small in (37r/7, 97t/7) . However, none of the transcendental behavior in 
(5.25) or (5.26) match with multiples of g\ and <72 directly. This is because there is an 
intermediate scale when |l + iqC| = Ca^ 7 ) where the effect of nonzero Ca enters 
into the WKB solutions (4.13) and (4.14) and the further approximation g\ and 02 as 
given by (4.17) are invalid. One can obtain the correct asymptotic behavior by going back 
to (4.13) and (4.14). It is easier in our case to directly obtain the same results from equa- 
tion (5.20) and (5.21). By linearizing about F\ = 0 , which is valid for large enough £1 , 
we obtain a pair of linear differential equations. The homogeneous part of those equations 
generate transcendentally small terms. By intoducing rescaled variables 

ft = /?i 9/7 ft (5.27) 
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into the resulting equations, they can be reduced to one homogeneous equation for the 
homogeneous solution F H (£ 3) . For £3 = 0(1) and >> 1 , this equation reduces 
to: 


d 3 F H 

~w 


l Pltl' 6 [1 + 0 ( 1 )] 


d 2 F H 

d £ 3 


- fl} e 3 3/2 [1 + 0(1)} 


dFjt 
d£ 3 


0 


The two independent WKB approximations to the independent solutions are: 


h. 


[/.*•«. [tV-Vt!'* + 36f s s/ ’] +•(!)] 


Fh 3 = e 


1 A 3 

7^1 


a, [tl^+y/tV + ^l' 2 ] + 0 ( 1 )] 


(5.28) 


(5.29) 

(5.30) 


For large £3 , the behavior in (5.29) and (5.30) matches with gi and <72 as given by 
(4.17) for small |l + i <?C| • On the otherhand, for small £3 , the behavior matches with 
F' H and F' H as given by (5.25) and (5.26) respectively. The o(l) in the expressions 
(5.29) and (5.30) are not important in our analysis other than that they are real for real 

6 - 

Including the leading order transcendentally small correction term, the asymptotic 
behavior of solution F\ for Arg £1 in the interval [0, 37t/7) is 

F 1 ~ — + •• + AiF I j l (5.31) 

On the other hand for Arg £1 in ( Zn/7 , 97t/7) , the transcendental correction cannot 
include a multiple of F Hl which gets large in this sector as £1 — » 00 and does not 
match with a multiple of (5.30) as it must in order to match with the behavior (5.2) in 
sector II of the ( plane in Tig. 3. 

As far as the differential equation (5.22), a unique solution is obtained by numerically 
integrating the solution from £1 — Le l3 *t 7 to £1 = £1 where L is chosen to be 12 
and £i = l with the condition (5.23) imposed as an initial condition. The choice of the 
initial point and initial condition was made to suppress F Hl where it is transcendentally 
large at £1 = £1 for large L . The solution was found to have a non zero imaginary part 
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at £1 implying that the coefficient of F Hl is not purely real implying that Im C\ in 
(5.1) is not zero. Thus no finger solution with smooth tip exists in this case. The results 
were checked by varying L and £1 . 

The leading order asymptotic solution Fi 0 is not valid for all £1 . Its derivative 
is singular at some £1 = £i 0 in the upper half £i 0 plane where the local singularity 
behavior can be directly deduced from (5.22) to be 

i 3/4 

K ~ -(4) (Cl-a „)- 3/4 (5.32) 

Thus the approximation F\ ~ F\ 0 becomes invalid near £ = £i 0 as the next order 
approximation F\^ becomes singular owing to inclusion of the ' term in (5.20). 
However, this does not affect the conclusion about the non-existence of solution. Here, we 
present the details of the analysis of this innermost region for the sake of completeness. 
We rescale variables once again in the neighborhood of £1 = £i 0 : 


(1 = £i 0 + /?r 12/7 & 


(5.33) 


Fi = ~/?I 3/7 F 2 (£ 2 ) 


(5 A) 


Then for £ 2 = 0(1) , to the leading order (5.20) and (5.^?) simplify to 


/ 


F 2 = 


+ -T7f 


n/ / 


F 1 V 3 F'J 


X 


X 


/ 


(5.35) 


jr S 

For large (2 > it is clear that the asymptotic behj^ r consisj^lt with (5.32) is given 

C x 

3/4 ^ 

F 2 = J£J (5.36) 


by 


The transcendental cor rection to *'hi§X av * or ^ oun( i by linearizing (5.35) about (5.36) and 
examining the homogeneous ’art/ e fi ua ^ on is a multiple of 


/ e a * (<) [1 + «(!)] 


(5.37) 


/ 
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This is transcendentally large for Arg £2 in (-4 tt/7, 4tt/ 7) and transcendentally small for 
Arg £2 in (4 tt/7, 127t/ 7) . It is necessary to suppress the transcendentally large behavior 
for large £2 when Arg & is in (- 47 r/ 7 , 47 t/ 7 ) since the transcendental behavior given 
in the above equation matches to the following transcendentally small correction for small 


/?! 


/‘ l Ui-rlJ ' 1 -Cl [1 +•(!)] 

3 'In 


(5.38) 


This was determined by linearizing ( r 20 ) and (5.21) about F\ 0 and looking for WKB type 
solutions for small f3\ to the resulting homogeneous equation. The expression in (5.38) 
matches to a constant multiple of (5.30) for large £1 and small £3 . On the otherhand, 
for large £3 and small |l + iqC, | , a multiple of (5.30) will match with C 2 g2 in sector 
2 in Fig. 3 as required. Since (5.35) is autonomous, the other degree of freedom of the 
solution Tn the second order differential equation is the arbitrariness of the choice of origin 
of £2 , but this does not~affect the matching. Thus a unique value of C 2 in (5.2) can be 
determined in this matching procedure. However, as we concluded before, the finger tip in 
this case is not smooth. 

»c. Case of Ca ^~ 1 (j_( 1-2A)~ 1 A 2 >> 1 or order unity. 

In this case, each of £2 , /?i are large, but /?3 can be large or small or 0(1) compared 
to unity. It is appropriat e f (- introduce t he transfor mation 


Fi(Zi) 


& = (h 3/4 r2/? (! + 

= -^ 3 / 2 r 4/7 C 1 + ? 2 ) 


2q< 


HO 


«i(*i) = - ap ?- 4/7 . (l + g } 


2 q 2 =iKL 


(5.39) 

(5.40a) 

(5.406) 


into equations (3.3) and (3.4), which to the ori^r in r 2/7 ? assumed small, 

for Cl = 0 ( 1 ) ou elimination of , redu ^ ^^^ 
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+ 


Pl l/ * 


FI" 


+ 3 (Ff - iy 


(5.41) 


l ( 4 i -- f , {) 3/2 2 ( ei -^ i ) 572 - 1 

Here, in this case we assume that Ca _1/2 A -1 (1 — 2A) >> 1 in order that 

r -2/7 ^3/4 >> i . Xhis assumption is found to be satisfied for the class of solutions 

found for Ca >> Q . This will be assumed to be the case. Note that for = 0(1) , 
Ca n x = O (A 1/2 (1 — 2A) _1/2 r 1/7 Ca fa 3/8 ) and using the definition of r and fa , this 
is found to 0(Ca 3/4 ) << 1 . Thus the Bretherton approximation used in the derivation 
of (5.41) from (3.3) and (3.4) is indeed valid. 

First, when Ca Q ~ 1/2 (1 — 2A) _1 A 2 >> 1 , it is clear that fafa^ 1 * << 1, 

Pi fa 7 ' ^ << 1 and 21/8 << 1. Of these small parameters /?i /? 2 7/4 is the largest. 


To the leading order 


F\ = F[ 


where F{ is determined by the nonlinear algebraic equation 


(5.42) 


K + ns = 0 


The next order correction is clearly 


F[ 


Pi 




11 pT 3(^i — F{ q ) 4/3 + 1 
Note that in the upper half plane, i r { i is singular when 

fa = 4 e i 3 "/ 4 3- 3 / 4 = fa 0 


(5.43) 


(5.44) 


(5.45) 


where the local behavior of ,F{ o and F{ i are given by 

f _ " e. t. \ t c~ c. \ 1 ! 2 


3 3 1/ 2 

- (g^io) (fa - fa 0 ) L/ ‘ 


n x 


pi i 


(5.46) 

(5.47) 


— o if 2 

as can be seen by analy zing (5.43) aftd (5.44). The minus sign preceding the ( ) term 

1/3 

in (5.46) is the only choice consistent with the behavior F{ o ~ —fa for large on 
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the real positive £1 axis (with the principal branch choice) as was determined numerically 
by solving the nonlinear algebraic equation (5.43) for F ^ for each £ , taking care that 
there is no large jump in the the argument of £1 — F{ g between neighboring £ points 
(assuring us that we are on the same branch of the Riemann sheet). Now, it is possible to 
continue to try to find the next order correction F\ 3 in a similar fashion. However, for 
our purpose it will suffice to note that for Ca >> Q (1 — 2A) X ~ ^ , the two term 
leading order behavior of F\ is: 

Fi ~ F lo + F i, = F a (5.48) 


To find transcendentally small corrections to this behavior, we linearize equation (5.41) 
about F a and look at the homogeneous part of the equation whose solution will generate 
terms with transcendental dependence in , /?2 and /?3 . By a standard dominant 
balance argument on the resulting linear equation, tedious, but routine calculation shows 
two possible behaviors for the solution to the homogeneous equation. These will be denoted 
by F Hl and F H , and their derivatives are given by 






4/3 F i 0 )‘ /S + l]+°( 1 )] + / 0 ’ [,«,-%■ >«/’ + . 

= tii - F u) e 


pi _ +.(i)] 

r H 3 ~ e 


(5.49) 

(5.50) 


Note that in (5.49), the occurrence of the second term in the exponent is due to the 
inclusion of F^ as given by (5.47). The product of and o(l) term appearing in 

the exponent, in (5.49) depends on the coefficients j3\ , /?2 and /?3 and terms of this 
product can in some cases be larger than unity and therefore exceed the second term in 
the exponent in (5.49) which_does not contain any parameter. However, it is convenient 
in our case to isolate the parameterlndependent 0(1) term in the exponent in (5.49), as 
we have done, and write the remainder which has no p aram eterjndependent 0(1) term 
as the product of and o(l) . For our purpose in determining^ the leading ordep— 

. ..... : jg g fe r - ' 
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selection rule, it will not be necessary to find any of the o(l) terms in the exponents in 
(5.49) and (5.50) aside from noting that they are real on the real £1 axis. To connect the 
behavior the behavior of transcendental terms in (5.49) and (5.50) to the far field behavior 
where transcendental terms are linear combinations of g\ and g*i as defined in (4.17), 
one needs to introduce an intermediate scale variable: 


£3 = /?r 12/7 i = /?r 12 / 7 /?f /4 ei 


= F B {i 2 ) 

Then for £3 = 0(1) , (3\ — * 00 , one finds the equation for F H reduces to 

d 3 h 


dti 


f - J tfer [1 + o(D I - Pi d ' 2 [1 + 0(1) 1^7 = ° 


dii 


The two possible solutons to this equation for large are 


F' Hl = e ' 


n , = « 


1 A 3 


/ 0 fS ^3 {t 1 ,'* W tl'* + 36 f’ /a ] + 0(1)] 


(5.51) 

(5.52) 


(5.53) 


] + «(!)] 

(5.54) 

] + °(1)] 

(5.55) 


For large £3 , the behavior in (5.54) and (5.55) matches with and g '2 respectively. 
On the other hand for small £3 , it is clear that (5.54) and (5.55) correspond to F' Hi and 
F' Hi as given by (5.49) and (5.50) respectively. The details of the Stokes lines in different 
nested inner regions £1 = 0(1) and £3 = 0(1) are shown in Figs. 5 and 6. The solid 
lines are the Stokes lines with respect to F ^ and the dotted lines are the ones respect 
to F' Hi . Thus, in order that there be no transcendentally large term in sector I and II 
of Fig. 3, it is necessary that the form of transcendental of transcendental correction to 
(5.48) for large £1 be given by a multiple of F' g in sector I and of F g ^ in sector II of 
Fig. 5. Also, note that g\ is completely real for £ on the imaginary axis between iq and 
i/q . Also, F' H ^ is real on the real and positive £3 axis. Thus the requirement that the 
imaginary part of the transcendental correction for / be completely real on some interval 
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on the imaginary £ axis including i , to the leading order, is equivalent to requiring that 
the multiple of g\ in the transcendental correction in sector I of Fig. 3 be real, which in 
turn implies that the coefficient of F^ in sector I of Fig. 5 be real. Now, to find what 
this coefficient must be, we introduce an inner variable in the neighborhood of £1 = £i 0 , 
where the approximation F' ~ F[ o becomes invalid. 


n = + 


■02 


01 

, 7/4 


1/3 


1 02 "* 7 X4 

Then for <,2 = 0(1) , to the leading order (5.41) reduces to: 


/ 1 \ 2/3 / 3\ 2/3 

(jft.) ( 2 ) F 2&) ( 5 - 5? ) 


+ F? = -F'{ + 


(5.58) 


where 

* = 3 (|) 1/3 (i h.y ut K* 1 (5.59) 

Under the conditions of this case, we note that /?4 << 1 . To the leading order, if we 
ignore the /?4 term completely, we get the exact solution 



Ai'(t 2 ) _ 

nm = 1 


(5.60) 


and this matches with (5.46) for |f 2 | >> 1 for argument £2 in the interval (— tt, tt/3) cor- 
responding to small |£i — £i 0 1 with Arg (£1 — £i 0 ) in the range (— 37 t/4, 77t/ 12) as can 
be seen from the asymptotic expansion of the Airy function and properties of its analytic 


continuation in different sectors over different ranges of arguments. To the behavior (5.60), 


one can add a transcendentally small correction for small /?4 or large £2 that is a multiple 
of e p * (ta-fa, ) awa y from the immediate neighborhood of the set of £ 2 , > the j th zero of 
At (£ 2 ) ? j ranging from 1 to 00 . For large £2 this correction is transcendentally small 
only when Arg is in the range (w/2, 37r/2) , i.e. Arg (£ i-£i 0 ) in (37r/8, 117 t/8) . 

This transcendental term matches with Fjj 3 . The multiplicative constant term is found 
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by requiring that the solution to (5.58) contain no transcendentally large term in £2 for 
Arg £2 in interval ( — 77 t/8, 7t/8) . However, we will not need this constant multiple in 
the determination of the finger width. The transcendental correction for large £ 2 with the 
argument of £2 in the interval (-7T, — 7r/3),i.e, Arg (£1 — fi 0 ) in (— 37 t/4, 7t/ 12) can 

be directly calculated from the behavior (5.60). We find that for Arg £2 in ( —7r > * 

F' „ - j ^' 2 e» £ * /a (5.61) 


and this matches with the behavior 


F[ = K + Al F 'hx 


(5.62) 


provided that 


Ai = 


. -1/2 [/ 4l ° + i] + “Ml d(l 

I llj 112 « 1 0 1 0 


r i *tt 


>‘ /S 


* o 


- T In (-^0) 


(5.62) 


where 


u 1 = 


U2 = 




ft 7 ' 4 ' 

ft ' 
■ft 7 ' 4 ’ 


4 v 2> 

1,3 (Knr 


On numerical evaluation, we found 

0 


Im j d(i Ui-i r i 0 ) 4/3 + 3 - _0 - 


.473 


Jm 


y ft 0 

/ 

0 


4(ei-fQ 1;3 ( 1 -^) L_ 

[3({i - FiJ 4 ' 3 + 1] 2(£l-£l.)J 2 


- ^ In (-fl.) 


(5.64) 

(5.65) 


(5.66) 

-1.27 

(5.67) 


The condition that A\ be real is then equivalent to 

O .473 1^2— [l + o(l)] - 1.27 + \ Arg «i = + Arg u 2 + n n (5.68) 

Pi 1 L 


35 



Thus, if we drop the o(l) term in (5.68), we obtain 


o /? 2 /4 {n 7T - | + 1.27) 
3 Ti 073 


(5.69) 


Where n is a positive integer. Then going back to the assumption made in this case, 

n must be a large integer. Generally, there is not much point in including the term 

-7t/ 8 + 1.27 in the right hand side of (5.68) because the product of /?J /4 //?i and o(l) in 

(5.68) that is neglected in (5.69) is larger than unity. However, we still include it to make 

the formula more accurate in the special case when this product is smaller than unity as 

is the case when each of the conditions 7^7 - << 1 and ^2 21/8 ^2 7/4 /^i « l,i.e. 

P 2 Pi 

for large but not too large n . Equation (5.69) implies that for jjrzrj frfc » 1 , for 
sufficiently small Ca , 


a 2 (i-a) 

(1 — 2A) 


C 1 ' 2 

~ 2.06 ( n + 0.279) - 


(5.70) 


The above expession is only valid for large n ,i.e. higher branches of solution since for 

small n , this would imply (1- 2a) e>/s — 0(1) • 

For the lower branches, i.e. when n is order unity, -^77 = 0(1) , implying 

1 

Ca G ~ 1 ' 2 (1 — 2A)~ 1 A 2 = 0(1) one has to resort to numerical calculations of the 

solution to the dilferential equation 


Fl + (a-F {) 1/3 

with the requirement that the behavior of 
( — 7r / 12, 7 tt/ 12) be given by 


fh 


1 — F'y 


(5.71) 


F{ for large £1 with Arg £1 in the interval 


Fi 


+ 


Pi 


(5.72) 


rf /3 

and that it be real on the positive real £1 axis for sufficiently large £1 . Such solutions were 
found numerically by integrating the first order equation (5.71) for F[ from Z e’ to 
£l = £l , where we chose L = 12 and £1 = 0 . with the asymptotic condition (5.72) 
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satisfied at the far end point. We required that Im F{(£i) = 0 in order to determine 

Sj— , We doubled L and changed £1 to other real positive values without noticing 
3 J 

any change in the final results to the accuracy quoted. We found the first two values of 

-£77 -r to be 9.848 and 15.69 implying 
3/s!' 4 


* 2 (1 ~ A) 

(1 — 2A) 


3.061 


g 1/2 
“c7 


(5.73) 


for the first branch compared to 


A 2 (l-A) 
(1 — 2 A) 


4.878 


gl/2 

Ca 


(5.74) 


for the second branch. If we were to base the prediction on the formula (5.70), then for 
n = 1 and 2, the coefficients on the right hand sides of (5.73) and (5.74) would have 
been 2.634 and 4.695. Though formula (5.70) is formally invalid for small n , for n = 2 , 
already we have a 4 percent accuracy of the asymptotic formula (5.70). 

If we compare with the experimental situation where the gap to width ratio e = 37 
and Ca are the two convenient control parameters, we have the prediction that for 


1 >> Ca » e, 


A 2 (l - A) 
(1 — 2A) 


= 2.776 



for the n = 1 branch and for n = 2 , 


(5.75) 


A 2 (l- A) 
(1 - 2A) 


4.424 



(5.76) 


Note that for 1 » Ca » £ 1/2 , i.e. 1 >> Ca » e 2/ 3 , from formulas (5.71), 
(5.73) and (5.74), for fixed integer n , A approaches zero like Ca i/* or Co s/4 • If on 
the other hand Q << Ca « C?^ 2 << 1 , i.e. for e << Ca << e 2 ^ 2 > A — ^ is 

• 1 C cl 

proportional to ^777 or — - — 

The most interesting case is when Ca = 0{ Q ^ 2 ) , i.e. Ca = e 2 ^ 3 because in 
that case the selected A is in the open interval (0, ^) and the different branches, i.e. 
different n , do not asymptote to the same limiting A . This means that if one of the 
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solutions were stable and the others unstable, there would be no mechanism for nonlinear 
stability as put forward by Bensimon. 

6. Case of A > 5 

For A > ^ , i.e. p 2 > 0 , we will restrict the analysis to the case of p << 1 , i.e. 
A close to one half. Our primary concern in this section will be to confirm the validity of 
the previously calculated selection rule (Combescot et al (1986), Shraiman (1986), Hong k 
Langer (1986), Tanveer (1987b), Dorsey k Martin (1987b) and Dorsey k Martin (1987)) 
for Ca « Q 7/ 3 . In this case, it is appropriate to introduce the following transformation 
of dependent and independent variables: 



(6.1) 

f = p 2 d(o 

(6.2) 

g = p 2 u,(0 

(6.3) 


Then to the leading order in p , equations (3.3) and (3.4) reduce to 


where 


f P' - 2£ 2 

(£ 2 - 1 -£D') 3/2 

IT / 3 




$2 


(s — 1 — £ D ') 


^176 = ~ iS ( D + 


^ du 

^ = ~ e 


S3 


(£ 2 -l-£D') 


a 1/3 


+ e l 


r/6 & 


{i 2 D" + ££>' - 2e 2 ) 


(| 2 -1 -ZD>) 


,ai/6 


r _ P 3 4 

G 7r 


62 = 3.878 


2 !0/3 ^1/ 3 


7 r 


v 12 


Q - 1 ' 2 Ca 1 ' 6 


S3 = 


1.3375Ca 2 ' 3 


2 2/3p8/3 

= !- 3375 ^3 6 1/2 Co 7 ' 6 ; 


,11/3 


(6.4) 

(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 


Note the differing definition of 6 from previous work (Tanveer, 1987b). This corresponds to 
multiplicative adjustment of the parameter G by ^ as would be necessary for the Mclean- 
Saffman theory to be valid even for Ca « G 3 (assuming zero contact angle between the 
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advancing interface and the gap plates) as noted earlier. The two equations (6.4) and (6.5) 
can be conveniently replaced by a single equation by defining a new dependent variable 


G(0 = £ 2 -l -££>' 


( 6 . 10 ) 


Then eliminating u between (6.4) and (6.5), we get 


G" + 


3 G' 2 
2 ~G 


- h £lle^L a ' = 


— i 6 64 


e ix/ 6 Ql 


( 6 . 11 ) 


( G 1 ' 3 

One needs to look for solutions to (6.11) for different possibilities of the parameters 62 , <53 , 
£4 and 6 such that they match with the outer solution in sectors I, II and III of Fig. 3. 
Depending on the relative size of these parameters, there is need for a set of nested inner 
regions as in section 5, In this paper, we do not address the possibility of solutions for 
different ranges of parameter. We only concentrate on the specific case Ca << and 

6 = 0(1) . It is clear that in this case 62 , 63 and £4 are each much smaller than unity 
so that the thin film effect totally drops out. In that case, in order to relate to earlier work 
(Tanveer, 1987b), it is convenient to go back to the original equation (6.4) and (6.5) which 
is now reduced to one equation 


Z 2 P" + ££>'- 2£ 2 
(£ 2 - 1 -tD'f 2 


-\6 D 


( 6 . 12 ) 


which is the same nonlinear equation as (Tanveer(1987b)) (equation 127 of that paper 
when — 6 and a = 1 , i.e. symmetry is assumed) as for the MS theory. From 
previous work, one needs to find solutions of (6.12) so that 


D ~ 



(6.13) 


for argument £ in the interval (— 7r,0) in order to match with the regular perturbation 
series in Q . For symmetric fingers, it is enough to require that the asymptotic behavior 
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(6.12) hold for large £ with Arg £ in the interval [— 7t/2, 0) with the requirement 
that the solution be completely real on the negative imaginary £ axis for sufficiently 
large £ . This condition determined 6 when it was order unity. For large <5 asymptotic 
analysis is possible and the details are given in the previous work (Tanveer, 1987b) which 
is equivalent to earlier analysis of the MS equations by Combescot et al. 

When 6 = 0(1) and Ca = 0(Q 7/3 ) , the 62 factor is order unity, whereas each 
of 63 and <$4 is negligible compared to unity. In that case, we can still reduce to one 
equation 


ep" + (D’ - ,. /6 


— e 


h 


= -i6 D 


(6.13) 


(£ 2 — 1 — £ D ') 3/2 " (£ 2 - 1 -£ D ') 1/6 

and numerical computation of these solutions using the same procedure as detailed in the 
earlier paper (Tanveer, 1987b) shows that for 62 = 0(1) , as we increase 62 from 0 to 
10, the corresponding 8 increases monotonically. This means that the finger is fatter then 
for the MS boundary conditions when parameters are such that the variation of transverse 
curvature is more important than the thin film leakage term. This is consistent with 
direct numerical computation of Schwartz & DeGregoria (1987) where inclusion of just 
the transverse variation of curvature term while neglecting flow into the thinfilm region 
resulted in fingers with width more than \ in the limit of small Q . 

We leave the detailed analysis for other ranges of Ca and Q for A > ^ for the 
future. 

7. Analytic theory for arbitrary capillary number: 


Eliminating g between (3.3) and (3.4) one finds 


W = Qh\ 


1-7 


4 QCa 


.0' 


{Ca n x ) n' x + £ h ' 2 


~ - [5(0 


k 1 (Ca 


i K V {Can x ) - / 13 (C) + 


7 r R 


^m°(Ca n x ) — /14(C) 


+ 


C 2 -l 

I \ 2 QCa 1 (h — h\ i (r , A 

V ~~^2 r W2 _ 1 ) m ( Ca n *) + a 


(7.1) 
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where n x , , l\ and I 2 are as in (2.18), (2.19), (2.20) and (2.21) respectively; hi , 
h2 , h 3 and /14 are as defined in (3.5)-(3.8). Note that on taking the derivative of (2.18) 
with respect to £ 


n, = — 


lo 1/2 Vo 


l 2 l l 
. 3/2 
‘1 


2 L ,1/2 

and this involves the second derivative of / 
For <5 = 0, 


j- 1/2 1 / 
<1 

.1/2 

*2 


+ 


, 1/2 » 

l l l 2 

. 3/2 


(7.2) 


C fo + MO - a 0 - m°(Ca = 0 


(7.3) 


as can be obtained from (7.1) or through direct analytic continuation of (2.23) where 


H(C)=/(mV« 2 C 2 -! 0 ) 

(7.4) 

i 2 . = i-pV + c(c 2 -i )/6 

(7.5) 

'1. = C 2 - p 2 + i (1 - C 2 )/«i) 

(7.6) 


Equation (7.3) has no exact solution. We can determine fo for |(| < 1 using the 

convenient power series representation in (2.24). However, we are interested in the behavior 
of fo and the next order regular perturbation term fi for |£| > 1 to find possible 
sources of nonuniformity of the expansion. If we define 


/i(c) = m- 1 ) 


(7.7) 


then numerical calculation of the coefficients in (2.24) allows us to calculate f' K and its 
derivative (which is related to /q(1/£) ) conveniently. Similarly the power series in (2.23) 
allows us to calculate h 3o . Equation (7.3) then becomes a nonlinear algebraic equation to 
determine /6(C) for |C| > 1 , which can be determined numerically by Newton iteration. 

Note that this kind of trick can be of wide applicability in other problems such as 
Kelvin-Helmholtz and Rayleigh- Taylor problems where use of numerical calculations and 
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conversion of nonlocal equations into local algebraic and differential equations can be 
exploited to study the evolution of singularities in the unphysical plane using a numerical 
procedure. 

Ass um ing f' K , /" and hz 0 as given, the nonlinear algebraic equation (7.3) 

determining /(, for any given C outside the unit circle, can be symbolically expressed as 

E(f&, 0 = 0 (7.8) 

where £ is a function of two independent variables. From implicit function theorem, 
(7.8) determines /6(C) such that C), C) = 0 . 

From equation (7.1), one can easily see that the next order perturbation term f\ as 
expressed in (2.25) satisfies an equation that can be written symbolically as 

E n (/(.(C), 0 /i(C) = fcl + n 2 fS (7.9) 

The expressions for 7Z\ and %2 can be obtained from (7.3) by a straightforward though 
lengthy algebra. Note that each of R\ and f?2 involve /6(C) beddes f' K f R , 
/1(C _1 ) and hz B . Each of the latter four functions are considered known as they are 
readily determined by using the power series representations for /o and f\ as |C~ 1 | < 1 • 
It is clear now that there is a singularity of /{ at C = Co where 


E n (/6(Co), Co) = 0 

(7.10) 

It is not difficult to see that at C — Co , where Ey' (/o(C), C) 

= 0 , the second derivative 

/6' is also singular. If we define 


C5 

ll 

(7.11) 

Ai = Rz(/o(Co), Co) 

(7.12) 

= E ( (/J(Co), Co) 

(7.13) 

■43 = E,.J. (/J(Co), Co) 

(7.14) 
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Then it is not difficult to see that in the neighborhood of ( = (o 

fi ~ A - Jt 3 (c - Co>1 ' 2 <715) 

Near £ = £o > which is a point of nonuniformity of the outer perturbation expansion, 
it is appropriate to introduce inner variables 


C - Co = «i f 


/' = mo) + u 2 F'(0 


where 


and 


u i = 


Af 6 G 1 ' 3 

a 1 J z aV z 


u 2 


a } 13 aV 3 s 1 ' 6 


aI /3 


Then for | = 0(1) , equation (7.1) to the leading order in Q reduces 


-F' 2 + £ = F" 


(7.16) 

(7.17) 


(7.18) 


(7.19) 


(7.20) 


which is exactly the same equation as in section 5c for £2 = 0(1) . Note that the highest 


derivative term F'" 
solution to (7.20) is 


does not appear at the leading order. However the appropriate 


„ _ Mf) 

MO 


(7.21) 


which is singular at the zeroes of Ai(£) and there exists an inner neighborhood around 
each of the zeroes that where the third derivative term will be important. 

Note in order that the behavior of F' be given by the specific solution (7.21) to 
(7.20), the structure of the Stokes lines should be qualitatively the same as that found in 
case (5c). This is expected to be true when Ca is less than 0.2 or so, i.e. within the 
experimental range; however this needs to be verified numerically. 
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The solution (7.21) is matched to the outer solution 


r ~ fb + S 1 ' 2 f{ + ..Cl (7.22) 

for large £ with corresponding £ on the imaginary axis and by following the procedure 
in section 5c, it is not difficult to see that the condition of C\ being real is equivalent to 

-g- 1/2 Im [WoiCi) - ^o(Co) + Q 1/2 m(Ci) - Wo(Co))] - 7T/2 + Arg U 2 

1 

- - Arg ui = n n (7.23) 

where £/ is any point on the imaginary £ axis just beyond £ = i , and n is some 
integer. 

Equation (7.23) is the selection rule. In order to get concrete results one needs m° , 
m 1 , ac° and /s 1 for complex arguments since the analysis in section (5c) suggests that 
£o will not be on the imaginary £ axis, where n x is real, but somewhat off the axis. 
From calculation, we find 

E n (/6(C), C) = C (1 + Q\) (7.24) 

where 

Q\ = - m°(Ca nj (()) - Ca - ( ! 1 - 'l.) ™°' (Co n“) (7.25) 

l 2 0 h 0 

Thus, the singularity point is where Q$ 0 = — 1 . For small Ca , noting the behavior of the 
Reinelt functions, this can only happen close to where Z^ 1 is large, i.e. near £ = i , the 
zero of the derivative of the conformal map corresponding to the Saffman- Taylor solution. 

The effect of finite Ca would perhaps be to move £o even further from the imaginary 
axis. However, there is an important issue that we are unable to answer at this point. If 
£o moves way off the imaginary £ axis, the determination of finger width selection involves 
knowledge of the analytic functions m° etc for complex arguments whose imaginary 
part can be large. Numerically, one only finds approximate solutions for m° ,.. for real 
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values of the arguments. There could two analytic functions that are quite close to each 
other on the real axis but differ significantly as we move off the real axis. If the problem is 
well posed for small Q , the deviations of the two functions just mentioned in the complex 
plane should not matter. However, if Co has a significantly large real part, it will matter. 
Does this mean that Co cannot move too far from the imaginary axis? Or does it mean 
that the problem is ill posed for even non zero Q « 1 for Ca = 0(1) in the sense 
that insignificant changes in the form of the functions , k ® and k* on parts of 

the real physical domain cause large changes in the selection. One is left to wonder if this 
has anything to do with the evolution problem in the time dependent case where fractal 
like structures (Maxworthy, 1987) have been observed in experiments at sufficiently large 
Ca . 


8. Summary of theoretical results and comparison with experimental data: 

In this paper, we have found concrete results only in the limit where both Q and 
Ca are small, i.e. e 2 << Ca « 1 , where e is the gap to width ratio £ . 

For A > ^ , we have found that the finger solutions of the MS theory where 

A ~ \ + constant Q 2/3 persist for Ca << Q 7/3 << 1 provided we make a 

multiplicative adjustment of the previously calculated constant (Tanveer, 1987b, Dorsey 
& Martin, 1987) by a factor of (|-) 2/3 . The possibility of other solutions with A > ^ for 
other ranges in the parameter space has not been investigated. 

For A < \ , for Q « Ca « 1 , we predict a discrete set of finger solutions for 
which 


where 


and for n >> 1 , 


A 2 (1 - A) Q 1 ' 2 

(8.1) 

(1 — 2A) “ Ca 

fcl = 3.061 

(8.2) 

k 2 = 4.878 

(8.3) 

k n ~ 2.06 (n + 0.279) 

(8.3) 
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Unfortunately, for most of the experiments to date, strict quantitative test for the 
theoretical result (8,1) is not possible because most of the experimental data shows that 
e is not small enough to make Ca and Q small at the same time (note Q — e 2 n 2 / (12C a)) . 
However, for the Tabeling-Zocchi-Libchaber experiment, the data shows general qualitative 
agreement with (8.1) with n = 1 , i.e. the first branch, as shown in table 1 for some of the 
smallest values of Q available. The theoretically predicted A in table I was obtained by 
numerically solving the nonlinear equation (8.1) using standard Newton iteration. For the 
range of the Tabeling et al experiment, we do not have any data for which the right hand 
side of (8.1) is small compared to unity. Note that for small but fixed values of Q as e is 
decreased, both the theory and experiment suggest that A approaches 1/2. From (8.1), 
we get 1 — 2A proportional to e 2 /(/ 3/2 for Q « e << Q 3 ^ 4 << 1. 

Table 2 shows that there is very little agreement with the Kopf-Sill & Homsy exper- 
iment where very skinny fingers were observed. It appears that in their experiment, the 
SPHR conditions will be invalid in some small neighborhood near the tip, because the 
product of the gap width b and the tip curvature (as calculated from the Saffman-Taylor 
theoretical formula) is order unity. It is clear that in the Reinelt(1987a) derivation of the 
interfacial boundary conditions, one assumes that the product of lateral curvature and 
the gap width is small everywhere. However, this product at points away from the tip 
is small for their experimentally observed finger and so the SPHR conditions hold except 
right near the tip. We suggest that this in some sense is equivalent to a tip perturbation 
on the regular Saffman-Taylor finger and that might explain the similarity of the observed 
features with that of Couder et al (1986) where a small bubble near the tip of the finger 
is found to dramatically affect the selection mechanism. However, we are not sure how 
to explain the observation of Kopf-Sill & Homsy that the skinny fingers are only found 
for extremely clean systems. The theory presented here assumes the SPHR conditions to 
hold everywhere on the interface and is therefore unable to account for the experiment. 
Note that in our theoretical prediction, when C? 1/2 << Ca << 1 , the width A is 
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small. However, the product of tip curvature and the gap width, which is proportional to 
A -2 e for small A (using Saffman- Taylor formula), is proportional to Ca 3 ^ 2 and is far 
smaller than unity as required for the validity of the SPHR conditions everywhere on the 
finger boundary. 

We conclude this section by discussing a finding of some theoretical importance. If we 
go back to the analysis of section 5c, we notice that it is necessary to carry out matching 
in a sequence of nested inner regions £3 = 0 ( 1 ), £1 = 0 ( 1 ) and then £2 ~ 0 ( 1 ). 
Note that the leading order equations in each of the inner regions contain parameters. If 
at the outset, we had linearized the equations about F = 0 , these set of nested regions 
would not exist. Indeed, it is the nonlinearity of F\ 0 that generated a nonuniformity 
in F\ x , which accounted for the inner region where £2 = 0(1) . The matching of the 
behavior of this region to the next outer region where £1 = 0 ( 1 ) is what determined 
the selection rule ( 8 . 1 ). A linearized equation would have failed to predict even the correct 
scaling. It may be pointed out that the linearized analysis as of Shraiman (1986) has 
been extensively used in selection problems both for the Salfman-Taylor problem as well 
as in dendritic crystal growth. For the Salfman-Taylor problem with the MS boundary 
conditions, the results are quite close to what one finds from a fully nonlinear analysis 
(Combescot et al (1986), Tanveer (1987b), Dorsey & Martin (1987)), the only difference 
being a small error in values of scaling constants. The reason that the linearized analysis 
works for the MS boundary conditions is that in the inner region, where the nonlinear 
analysis is relevant, there is only one parameter in the nonlinear equation that 

has to be determined by matching to the outer solution. If one replaced this nonlinear 
equation by a linear one, there will still be the same one parameter in the problem though 
the actual numerical value of that parameter will not be correct. However, this parameter 
contains all the scaling information. This is obviously not true in our problem with the 
thin film and the nonlinear analysis is absolutely essential to get the correct scaling laws 
as well as numerical constants. 
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We have presented an analytic theory for the selection of Saffman- Taylor finger in the 
presence of the thin film that is neglected in the previous theory based on the Mclean- 
Saffman boundary conditions. Precise scaling laws are calculated. There is prediction of 
very skinny fingers. Unfortunately, direct quantitative verification from existing experi- 
ments is not possible because Bretherton’s results are only valid for rather small Ca , 
and therefore a very small gap width to cell width ratio, far smaller than the existing 
experiments, is necessary in order that Q be small enough for the validity of our scaling 
laws. However, there is qualitative agreement with one set of experiments. The completion 
of details for arbitrary Ca provides some exciting possibilities that will be left for the 
future. The importance of nonlinear analysis has been pointed out. The linear stabil- 
ity of the steady states for boundary conditions incorporating thin film effects is also an 
important problem that is left for the future. 
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e 

Ca 

l/B 

G 

3 061 ££ 

Expt. A 

Theor. A 

1/112.5 

0.0395 

6000 

0.165 x 10" 2 

3.14 

0.48 

0.481 

1/65 

0.118 

6000 

0.165 x 10" 2 

1.05 

0.47 

0.447 

1/112.5 

0.0263 

4000 

0.247 x 10" 2 

5.77 

0.49 

0.489 

1/65 

0.0789 

4000 

0.247 x 10' 2 

1.93 

0.48 

0.470 


Table 1. Comparison with Tabeling et al (1987) experiment for a few values of e and 
Ca for which Q was small. Note that B = G/ir 2 , and is the same B introduced by 
Tabeling et al. Also, note that values of experimental A quoted here are eyeball readings 
from Fig. 8 of the Tableing et al (1987) paper. 


e 

Ca 

G 

3.061 

Expt. A 

Theor. A 

j 0.38 

X 

10- 2 

0.25 

0.49 x 10“ 4 

0.089 

0.15 

0.24 

O 

CO 

oo 

X 

io- s 

0.054 

0.22 x 10" 3 

0.85 

0.20 

0.44 

0.38 

X 

lO" 2 

0.0128 

0.95 x 10" s 

7.35 

0.25 

0.49 

0.38 

X 

To^" 

0.0036 

0.33 x 10“ 2 

48.44 

0.31 

0.50 


Table 2. Comparison with Kopf-Sill k Homsy (1987) experiment. 
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Figure 2. 



54 


Figure 3 
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Figure 6. 
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